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Abstract 

Paired states, trions and quarteting states in one-dimensional SU{4) attractive fermions are 
investigated via exact Bethe ansatz calculations. In particular, quantum phase transitions are 
identified and calculated from the quarteting phase into normal Fermi liquid, trionic states and 
spin-2 paired states which belong to the universality class of linear field-dependent magnetization 
in the vicinity of critical points. Moreover, unified exact results for the ground state energy, 
chemical potentials and complete phase diagrams for isospin S = 1/2, 1, 3/2 attractive fermions 
with external fields are presented. Also identified are the magnetization plateaux of = M^/S 
and = 2Ms/3, where Mg is the magnetization saturation value. The universality of finite-size 
corrections and collective dispersion relations provides a further test ground for low energy effective 
field theory. 

PACS numbers: 03.75.Ss, 03.75.Hh, 02.30.IK, 05.30.Fk 
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Experimental advances with higher spin fermionic systems of ultracold atoms present a 
unique opportunity to rigorously test the current understanding of molecular superfluids 
and more generally to probe the nature of quantum many-body systems. Three-component 
ultracold fermions give rise to a phase transition from three-body bound states called trions 
,n.o .Ke BCS pa,„,. ..e Q ^ I. pa.,cu,a, .,„-3/2 ,n..ac«n. a.^.c fe...on. have 

nnn 

been predicted to exhibit a quarteting phase, i.e., a bound state of two BCS pairs [2, 14, 15[, 



and BCS pairing with total spin-2 

considerable interest in the context of the one- dimensional 
ultracold atoms 



-ICOL Hi U 



a]. Such molecular superfluids have currently received 

ID) lattice Hubbard model of 



11| due to new progress in 



and degenerate quantum gases j^. 
experimental realization of highly degenerate atomic gases |12( |. 

Furthermore, recent experiments on systems of ultracold atoms confined to one dimension 
(ID) [l3, 14, 15] have revived interest in Bethe ansatz (BA) integrable models of interacting 



bosons and multi-component fermions. The nature of "no diffraction" in the many-body 
scattering matrix of ID integrable models results in key features of quantum many-body 
physics which are specified by the phenomena of spin-charge separatioii, scaling dimensions 
and universality classes of quantum phase transitions and criticality 18] . In this 

Letter, we study complete phase diagrams and quantum phase transitions in integrable ID 
SU{A) attractive fermions with external magnetic fields. We demonstrate that quantum 
phase transitions from quarteting states into phases of normal Fermi liquid, trions and spin 
paired states are fully controlled by Zeeman splittings. In particular, the unified exact results 
obtained for isospin S = 1/2, 1, 3/2 attractive fermions display a simplicity and universality 
which gives insight in understanding high spin paired states and spin liquid behaviour in 
multi-component interacting fermions. 

The model. We consider a 5-function (contact potential) interacting system of 
atomic fermions with equal mass m which may occupy four possible hyperfine levels 
i = 1, ■ ■ ■ ,4) labeled by isospin states ±3/2, ±1/2 and, constrained by periodic boundary 
conditions to a line of length L. The Hamiltonian 20 
kinetic energy Ho = -^Y^iLi -£2, interaction energy Hi = gmYl 
Zeeman energy = Ylt=i Here A^* is the number of fermions in state \i) with 

Zeeman energy (acting as a chemical potential /li [3l]) determined by the magnetic mo- 
ments /i^ and the magnetic field B. The spin-independent contact interaction qid remains 
between fermions with different hyperfine states and preserves the spins in each hyperfine 



2l| reads H = Ho + ni + E, with 
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states, i.e., A^* with i = 1,...,4 are good quantum numbers. Although these conditions 
appear somewhat restrictive, utihzing the broad Feshbach resonances, it is possible to tune 
scattering lengths between different states close to each other to form high degeneracy Fermi 



gases [12 



13|, 



22l |. Thus the model still captures the essential physics relevant to the multi- 
ple phases of molecular superfluids. The coupling constant gu) = —h^c/m with interaction 
strength c = — 2/aiD determined by the effective ID scattering length oid 23||. For simplic- 
ity, we choose the dimensionless units of ^ = 2m = 1 and use the dimensionless coupling 
constant 7 = c/n with linear density n = N/L. 

For an irreducible representation [^4^43^32^2 j^A^ij ^ four-column Young tableau encodes 
the numbers of unpaired fermions (A^^i), bound pairs (A^2); trions (N^) and quarteting states 
(A^4) with Ni = — N^~^^ and A^^ = 0. For convenience in calculation, we rewrite the 
Zeeman energy as Ez = — Yl^=o HiNi. Here Hq = J2t=i ^1/4 is irrelevant because Nq = N 
is fixed. The other values Hi,H2, are related to Zeeman splittings Aj+i j 



with 



i = 1,2,3 via the relation 







f' 


2 


i\ 


/a 




1 










H2 


~ 4 


2 


4 


2 


A 






V 


2 




u 



^21 

U3 



(1) 



We shall find that equally spaced (linear) Zeeman splitting, i.e., Aa+ia = A for a = 1,2,3, 
drive a smooth phase transition from a quarteting phase into a normal Fermi liquid, as 
depicted in part (A) in Fig [H Unequally spaced (nonlinear) Zeeman splittings may trigger 
spin-neutral bound states which are illustrated in parts (B) and (C) in Fig [H It follows 
from the relation ([1]) that for linear Zeeman splitting, the magnetic moments of a trion, a 
bound pair and the unpaired fermion are §, 2 and |, respectively. The quarteting state 



remains a spin singlet. 

The energy eigenspectrum is given in terms of the quasimomenta kj of the fermions via 
E = Ylf=i which satisfy the BA equations 

A/i 

exp(iA;jL) 
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Mi. The parameters Aa"* with 
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^^-^^iic' 
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1,2,3 are the spin 



rapidities, where we denote A^'' = ka and Aa^ = 0. The quantum numbers are given by 

Charge bound states. For attractive interaction, the BA equations ([2]) admit charge bound 
states and spin strings. In particular, the SU{A) symmetry acquires three kinds of charge 
bound states: quarteting states, trions and bound pairs. The patterns of these bound states 
and spin strings determined by ([2]) underpin the nature of quantum statistics and many-body 
effects in the atomic system. 

In the weak coupling regime, i.e., L\c\ <^ 1, we find that the imaginary parts iy of the 
charge bound states are the roots of Hermite polynomials Hk of degree k. Specifically, 
Hk{^^y) = 0, with k = 2,3,4: for a bound pair, a trion and a quarteting state, respectively. 
This result is indicative of a universal signature of many-body cooperative effects driven by 
dynamical interaction. The real parts of the quasimomenta deviate slightly from the values 
determined by Fermi statistics for the c = case. In this weak coupling limit, the BA 
equations ([2]) reduce to Gaudin model-like BA equations [2J] which describe the multiple 
charge bound state scattering. Using these root patterns, we explicitly obtain the energy 
(in units of h?/2m) 

2 ^ 1 K K 

k=i j=i k=j+i i=j 

4=1 j=i+l j=l J 

This result unifies the ground state energy results for two-, three- and four-component 
fermions (where k = 2,3, 4, respectively) for weakly attractive interaction. In this result 
the densities of unpaired fermions, BCS pairs, trions, and qarteting states are denoted by 
na = Na/L with a = 1,2,3,4, respectively. The ground state energy ([3]) is dominated by 
the kinetic energy of composite particles and unpaired fermions and has a mean field theory 
configuration, where the interaction energy accounts for density density interaction between 
charge bound states and between charge bound states and unpaired fermions. For the weak 
coupling limit, spin- neutral bound states are unstable against thermal and spin fluctuations. 
They form a gapless superconducting phase. 

On the other hand, for the strongly attractive regime L\c\ 3> 1, the imaginary parts of 
the bound states become equal-spaced, i.e., a qarteting state has the form ki = A^^^^ ±i3|c|/2, 
A(^)±i|c|/2; the trion state is kj = Aj^^±i|c|, A^-^'* and for the bound pair kr = Ar"'^^±i|c|/2. The 
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corresponding binding energies are given by = c^i{i'^ — 1)/12 with £=1,2, 3, respectively. 
Substituting these bound state root patterns into the BA equations ([2]), we exphcitly obtain 
their real parts A^^^, with £=1,2, 3. This leads to a unified expresion 



+ + (4) 

k=l \ I I / 

for the ground state energy for two-, three- and four- component attractive fermions (in 
units of h^/2m) where the functions Ai = 4^2 + 2^3 + ^2 = 2ni + ^2 + ^ + 
A3 = ^ + ifa + n3 + 5||iandA4 = f + ^ + ^ + The functions reveal the 
scattering signature in different channels, for example, rii does not appear in Ai due to the 
lack of s-wave scattering for unpaired fermions. We note that for two-component attractive 



fermions the terms involving ris and n4 should be excluded [26|] whereas for three-component 
attractive fermions does not appear. The unified structure of the ground state energy 
dl]) can be amended with appropriate Af^ functions for higher spin fermions. 

Charge bound states in equilibrium. The BA equations ([2]) in principle give the complete 
quantum states of the model. However, at finite temperatures, the true physical states 
become degenerate. In the thermodynamic limit, L, N — »^ oo with N/L fixed, the grand 



partition function is given byZ = tr(e-^/^)=e-^/'^[l6| , [l9|, 1251] , where the Gibbs free energy 
G = E + Ez — liN — TS, the chemical potential /i, the Zeeman energy Ez = — X]i=o ^iNi 
and the entropy S are given in terms of densities of charge bound states and spin strings 
subject to the BA equations ([2]). The equilibrium states are determined by the minimization 
of the Gibbs free energy, which gives rise to a set of coupled nonlinear integral equations 
— the thermodynamic Bethe ansatz (TBA) equations. Following the TBA treatment for 



spin-1/2 attractive fermions [16(], we obtain the dressed energy equations 

-K,3,5]*e(3)(r)-[a2,4,6]*e(^Hr), 



e 



(3) 



(A) = 3(A2 - /i^) - a2 * e(^)(A) - [a,^s] * e^'\X) 



-[a2,4]*e(=')(A)-[ai,3,5]*e(^)(A), 
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= e -^^ -J2a,*e^'-^'\k), (5) 



e 



i=l 



which will be used to study quantum phase transitions at zero temperature. In these 
equations the function aj{x) = (^jcf-zf+x^ ^'^'^ * ^^"H^) = Ig^^ji^ ~ y)^^"'Hy)dy is 
the convolution. Furthermore, we have used the abbreviation aij^k = + clj + CLk- We 
denote the dressed energies e^'*-' with a = 1,...,4 and the effective chemical potentials 
fi^ = Ha/a + fi + c^{a? — 1)/(12) with = for unpaired fermions, pairs, trions and 
quarteting states, respectively. The negative part of the dressed energies e*^"-* (x) for x < Qa 
corresponds to the occupied states in the Fermi seas with the positive part of e^""^ corre- 
sponding to the unoccupied states. The integration boundaries Qm characterize the "Fermi 
surfaces" at e('")(Q„) = 0. 

The Gibbs free energy per unit length at zero temperature is given by G = 
X]m=i ^ I-Q ^^"^\x)dx. The dressed energy energy equations (jS]) describe the band fillings 
with respect to Zeeman splittings Hi and chemical potentials and provide complete phase 
diagrams and quantum phase transitions for the model. Solving the equations ([5]) by itera- 
tion among the relations — |^ = n, = Ui, i = 1,2,3 and the Fermi point conditions 
e^^^Qm) = gives the effective chemical potentials 

TT^ ^ \ \C\ J \C\ C2 ^ ^ 

which characterize Fermi surfaces of stable spin-neutral states and unpaired states. In 
this equation we have introduced an inner product ■ A with A = {Ai, A2, A^, A4) and 
Bh, - Here B4 - I — , ^, igg-j, B3 - I — ^j, B2 - 

if) ^^'^ ^1 ~ Ts)' '^^^ A-i are as given above and / is the identity 

vector. The expression ([6]) for the chemical potentials contains the previous results for 



isospin 5* = 1/2, 1 Fermi gases 0, 126I, l27|] . We have found that the energy (jl]) derived from 



the discrete BA equations for arbitrary population imbalances can also be obtained from 
E/L = fin + G + riiHi + n2H2 + n^H-^. This indicates that the BA spin-neutral states 
comprise the equilibrium stable states in the thermodynamic limit. 

Magnetism and quantum phase transitions. The low-energy excitations split into collec- 
tive excitations carrying charge and collective excitations carrying spin. This leads to the 
phenomenon of spin-charge separation. The charge excitations are described by sound modes 

riQ 

with a linear dispersion. The spin excitations are gapped [5|, 17|] with a dispersion relation 
^uip) = a/ + v^p'^ where Aj, is the excitation gap and v,^ is the spin velocity in spin branch 
u. However, for strong attractive interaction the low energy physics is dominated by charge 
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density fluctuations. This is because the spin wave fluctuations are fully suppressed by a large 
energy gap at low temperatures. This configuration is evidenced by the universality class 
of finite-size correction to the ground state energy E{L, N) — LE^ = — St=i "^k, where 
the central charge (7 = 1 for the U{1) symmetry and charge velocities Vk ~ '^^^(1 + 
with k = 1, . . . ,A are the charge velocities for unpaired fermions and charge bound states. 
For the singlet ground state, the spin velocity Vg ~ + sp^) the spin-3/2 hyperfine 

branch, which is divergent due to the energy gap. 

We find from the dressed energy equations ([5j) that quantum phase transitions driven 
by Zeeman splittings can be determined by three independent external field-energy transfer 
relations 

H^ = ^c^ + u' -u\ H2 = 2c^ + 2{u^ -u^), = '^c^ + 3{u^ - u^) . (7) 

The Fermi surfaces and charge bound states are fully controlled by the Zeeman splitting 
parameters. The complete phase diagrams are determined by the equations (I7j) and certain 
combinations of these equations. Without loss of generality, we consider only terms up to 
order of I/I7I in the following analysis. 

Using the energy transfer relations ([7]), we find that linear Zeeman splitting A lifts the 
SU{4) symmetry to f/(l)^ symmetry for I7I ^ 1, i.e., linear Zeeman splitting does not 
favor spin- neutral bound states (recall part (A) in Fig[T]). The lower critical field Ad ~ 
^ — ^^^(1 + isfr) diminishes the gap, thus the excitations becomes gapless. Using the 



definition of magnetization = M^/Ms with = |ni + 2n2 + fns and Mg = |n, we 



384 v-^ ' 18|7|' 

fwQ+^riTi — M^/A/T iinl-V, A/f2 — I I ^^^^ _ 

see that in the vicinity of A^i, the system exhibits a linear field-dependent magnetization of 
the form = 192(A — Aci)/{mT^) with a finite susceptibility y = dm^/dA ^ 192/ (rni^). 



13, 



18|. When the Zeeman 



This result provides a testing ground for low energy field theory 
splitting A is greater than the upper critical field Ac2 ~ ^ + ^"3""^ (1 — g^) the system is 
fully-polarized into a normal Fermi liquid. For the intermediate regime, A^i < A < Ac2, the 
quarteting state and unpaired fermions coexist. The phase transition at the critical point 
Ac2 belongs the same linear field- dependent universality class. A plot of the magnetization 
vs Zeeman splitting is given in Figure O 

For nonlinear Zeeman splitting the quarteting state can break into two spin-2 bound pairs 
as depicted in part (C) of Figure [H In order to trigger such a paired phase, we let A43 = A32 
in the relation ([1]). In Figure [3] we demonstrate the resulting interplay between the quantum 
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phases of quarteting states (phase Q), spm-2 pairs (phase D) and unpaired fermions (phase 
U). We see clearly that the quarteting states are stable for A21 < — -^^^(1 + ~ ^43 
and A43 < — ^§^(1 + slp^) " nf^- The phase diagram shown in Figure [3] is determined by 
the first two equations of (171) and A21/2 = c^/4+/i^— /i^ describing the mixed D+U phase. In 
order to clearly see magnetization plateaux, we choose a simple linear relation A21 = /A43. 
We then find that / < |(1 — j^) and magnetization plateaux = 0,2Ms/3, occur, 
where Mg = 3n/2 is the saturation magnetization (see Figure [2]). This is obvious because 
the lines fixed by the slope / simultaneously pass the Q, D and U phases. Moreover, if the 
Zeeman splitting A32 is large enough, the phases D and D + U form the same phase diagram 
for spin- 1/2 interacting fermions with polarization 27 1. 

If we set A21 = A32 in the relation ([T]), the Zeeman parameters may trigger a phase 
transition from quarteting states into trions (recall part (B) in Figured]). Figure H] shows 
an exact phase diagram in the A2i(A32) — A43 plane. Varying the Zeeman splitting A43 
and A21 reveals smooth phase transitions from quarteting states into trions (phase T) or a 
normal Fermi liquid (phase U). Similarly, we choose A43 = /A21 and see that magnetization 
plateaux = 0, Ms/3, Ms occur if / > |(1 + yI^)- The critical points for the plateaux 
can be analytically determined from the first and third equations in ([7]) and an additional 
equation A = + fi^ — fi^ for the mixed T + U phase. Actually, if A43 is large enough, 
the model reduces to attractive fermions with three hyperfine levels [9j. We point out that 
all phase transitions are of second order with a linear field-dependent magnetization in the 
vicinities of critical points. 

To conclude, we have presented unified exact results for complete phase diagrams and 
quantum phase transitions from quarteting states into spin-2 paired states and trions in 
ID SU{4) fermions with population imbalance. The ground state properties and magnetic 
effects provide a testing ground for low energy effective field theory and a bench mark for 
experiments with multicomponent ultracold fermionic atoms. 
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FIG. 1: Depiction of phase transitions from quarteting states into (A) normal Fermi liquid, (B) 
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FIG. 2: Magnetization as a function of the Zeeman splitting parameter A. The curves shown, 
obtained from ([7]) with the effective chemical potentials given in ([6]), are for linear Zeeman splitting 
(dashed line) with c = — 10 and n = 1 and unequally spaced Zeeman splitting (solid line) with 
c = —8 and n = 1. The 2/3 magnetization plateau occurs through subtle tuning of the Zeeman 
splitting parameters. Here, e.g., A21 = /A43 = /A32 with / = 0.36. 
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FIG. 3: Typical phase diagram in the A21 — A43(A32) plane with A43 = A32 in the strong coupling 
regime c = — 8 and n = 1. In this case we see clearly that Zeeman splitting can trigger a quantum 
phase transition from a quarteting state (phase Q) into two spin-2 bound pairs (phase D), as 
depicted schematically in part (C) of Figure [TJ Other phases involved are unpaired fermions (U) 
and mixed phases. 
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FIG. 4: Typical phase diagram in the A2i(A32) — A43 plane with A21 = A32 in the strong coupling 
regime c = — 5 and n = 1. In this case we see clearly that Zeeman splitting can trigger a quantum 
phase transition from quarteting states (phase Q) to trions (phase T), as depicted in part (B) of 
Figure [TJ Other phases involved are unpaired fermions (U) and mixed phases. 
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